S-wave TT — TT Scattering Lengths in the SU(2) NJL Model 

Beyond Mean-field Approximation 

Mei Huang^, Pengfei Zhuang^'^, Weiqin Chao^'^ 
^ Institute of High Energy Physics, Chinese Sciences Academy, Beijing 100039, China 

'g\'. „ 

Q> ■ China Center of Advanced Science and Technology, World Laboratory,Beijing 100080, China 

a^; 

^~^ • ^ Physics Department, Tsinghua University, Beijing 100084, China 

^ . (February 1, 2008) 

cn ; 

(N. 

. i The S-wave tt — tt scattering lengths ao and 02 are calculated to the lowest 

> : 

00 . order of l/Nc expansion in the general framework of SU(2) NJL model beyond 

^y" ', mean-field approximation. It is shown that using the universal curve of oq 

Qs I and 02 the four NJL parameters, i.e., the current quark mass ttiq, the four 

"^ I fermion coupling constant G, the quark momentum cut-off A^ and the meson 

I , momentum cut-off A^, also the S-wave tt — tt scattering lengths in NJL model 

Qh. 

r~| ! can be uniquely determined. 

>: 

X- 



I. INTRODUCTION 

Pion-pion scattering at threshold became more interesting recently, because it could 
provide a test to the mechanism of chiral dynamics. 

The Nambu-Jona-Lasinio (NJL) model []I| has been regarded as a cornerstone to under- 
stand the chiral dynamics by assuming that the spontaneous chiral symmetry breaking be 
triggered by the large quark condensate in the vacuum. And it is described through two 
Schwinger-Dyson (SD) equations for quark propagator and meson propagator respectively 
- 0] . The S-wave tt — tt scattering lengths oq and 02 have been investigated to the leading 
order of 1/A^c expansion (where Nc is the number of color degrees of freedom) in the mean- 
field approximation, i.e.,Hartree plus RPA, of SU(2) NJL model 0-0. We know that 
at the mean-field approximation level, the two SD equations are not fully coupled, i.e., the 
solution of the meson SD equation has no feedback to the quark propagator, which induces 



the loss of the information of meson modes ||T0|. Therefore, it is necessary to go beyond 
mean-field approximation. 

In this paper, we calculate the S-wave pion-pion scattering lengths a^ and 02 in the 
general framework of SU(2) NJL model including current quark mass explicitly [l^, which 
is based on a chirally symmetric self-consistent scheme considering meson cloud contributions 
01 . There are four parameters to be fixed, i.e., the current quark mass tuq, the four fermion 
coupling constant G, the quark momentum cut-off Aj and the meson momentum cut-off A^. 

It is very difficult to extract precise values of the two scattering amplitudes in the ex- 



perimental analysis |Tl|]. However, from either forward dispersion relations or from the Roy 



equations, one can find that ao and 02 are constrained to lie on a universal curve [|I2|. We 
try to find an appropriate series of NJL parameters by using the new universal scattering 
length relation newly evaluated by M. G. Olsson in ||13[ together with the two observables 
m^ = 139 MeV and U = 92.4 MeV. 



II. SU(2) NJL MODEL BEYOND MEAN-FIELD APPROXIMATION 

The two-flavor NJL model is defined through the Lagrangian density, 

L = ij{iYd, - mo)^ + Giii^tpy + {^ti,r^)% (1) 

where G is the effective coupling constant with dimension GeV~^, tuq the current quark 
mass, and ip,ip quark fields with fiavor, colour and spinor indices suppressed, assuming 
isospin degeneracy of the u and d quarks. 

We first briefiy review the general scheme of the NJL model beyond mean-field approxi- 
mation |T0|. Including current quark mass explicitly, the quark self-energy m expanded to 
0{1/Nc) order can be expressed as 

m = rriQ + rriH + Sm, (2) 

where mn and 6m are the leading 0(1) and subleading 0{1/Nc) contributions shown in Fig. 
1. The solid lines in Fig. 1 indicate quark propagator S{p) = l/(^ — m) with full m. The 
quark condensate < gg > is a one-loop quark integral 
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<qq>= — <i^i^ >= -^iN^m / j^^^^- ^. (3) 



The corresponding meson propagator DM{k) (M means vr or a) has the form 

-^^*^^^) = l-2GnM(A:)' 

UM{k) = nS^^)(fc) + 6ut]{k) + 6U^^]{k) + 6uf}{k), (4) 

where 11^ is the meson polarization function, which includes the leading order U\j and 
subleading order 511^/' , shown in Fig. 1. The above constituent quark mass and the meson 
propagator, namely the Eqs. (^) and (^, or the Feynman diagrams in Fig. 1, form a self- 
consistent description of the SU(2) NJL model to the subleading order of l/Nc expansion, 
which is different from the earlier calculations [|1 - [0 where only the mean-field quark mass 
rriH and the Random-Phase- Approximation (RPA) meson polarization function U^j are 
considered. 



The meson mass rriM satisfies the total meson propagator's pole condition 

1 - 2GUMik^ = ml) = 0, (5) 

and the meson-quark coupling constant gMqq is determined by the residue at the pole 

9Mqq = {dUM{k)/dk^)-^\f,2^^2^. (6) 

Another important quantity in the meson sector is the pion decay constant /^ which generally 

satisfies 

n on' ^ ' 

y-irqq ■^'-J 

In the chiral limit, fj^ satisfies the Goldberger-Treiman relation f^[k)g.,rqq{k) = m. 

III. TT-TT SCATTERING AT THRESHOLD AND THE UNIVERSAL CURVE 

Now we turn to the calculation of the S-wave tt — tt scattering length oq and 02. The 



invariant scattering amplitude can be generally written as |12 



Tab,cd = A{s, t, u)6abScd + B{s, t, u)5achd + C{s, t, u)5adhc, (8) 

where a, h and c, d are the isospin labels of the initial and final states respectively, and s, t 
and u are the Mandelstam variables, s = {pa + Pt)'^, t = {pa — PcY, u = {Pa — PdY, where 
Pa ~ Pd are initial and final isospin momentum. Using perfect crossing symmetry, one can 
project out isospin amplitudes 

Ao = 3A + B + C, Ai = B^C, A2 = B + C. (9) 

In the limit of scattering at threshold, ^/s = 2171^^, t=u=0, the S-wave scattering lengths a/ 
(given in units of m'^) are: 

aj = -^Aj{s = 4ml,t = 0,u = 0), 1 = 0,2. (10) 



It is known that the two scattering lengths ao and a2 at threshold are not independent, 
they are constrained to lie on a universal curve |T^ |T^ : 



12 f°° da 
2ao - 5a2 = — / — — — ^ImA/,=i(g), (11) 

vr Jo g(l + q^) 



where q is the cm. momentum in units of raj^. And in [|13[, M. G. Olsson evaluated the new 

universal relation as: 

4 4 

2(^0 - ^-Oq) - 5(a2 - --al) = Lq, (12) 

on OTT 

i'^-rTTT^l-^-l"). (13) 

vr Jo g(l + g"') 
The integrand of Lq = 0.58 ±0.015 m~^ in [|T^ includes the dominant contribution of p(770) 
and other small contributions of higher resonances. 

We will calculate the scattering lengths to the lowest order in l/N^ expansion of the 
process tttt -^ vrvr in the general framework of SU(2) NJL model beyond mean-field approxi- 
mation. The Feynman diagrams include box and a— exchange diagrams, which are the same 
as those in [0]. The only difference is that in 0] quark mass rriH is in the mean-field approx- 
imation and meson polarization function ^]^^ in RPA, while in our calculations, quark 
mass m and meson polarization function IIm include the subleading order contributions. 

To keep the lowest order of the a exchange diagrams, the internal a propagator should 
be in leading 0(1/A^c) order, which has the same form as that in RPA 11]^^ |T^. So 
what we need to do is using the expression formulae given in |^ directly, and replacing the 
quantities in the mean-field approximation with those beyond mean-field approximation, 
i.e., the external pion propagator becomes the total propagator Eq. (^), the pole for quark 
propagator becomes the total quark self-energy m defined in Eq. (^), and the coupling 
constant gT^gg now is the one expressed in Eq. (|). 

IV. NUMERICAL RESULTS 

For the numerical calculations, we adopt the external momentum expansion method 
as discussed in detail in |jlO|. We introduce a quark momentum cut-off Aj in Pauli-Villars 



regularization and a meson momentum cut-off A;, in covariant regularization for the divergent 
momentum integrals. Using only the two experimental observables m^r = 139 MeV and 
/^ = 92.4 MeV, one can not give fixed values of the four parameters in the model, namely 
the current quark mass tuq, coupling constant G, and the two momentum cuts Aj and A^. 
We introduce one more free parameter z = A^/Af, which characterizes the meson cloud 
contributions. Especially, in the limit of 2; = 0, the model goes back to the mean-field 
approximation automatically. 

For each z, with quark mass m changing from 200 MeV to 1200 MeV, we solve the three 
equations (0), (|) and (0) to get a series of A/, niQ and G, and then calculate other quantities 
like < qq > and g-^qq. Finally, we can calculate the S-wave vr — vr scattering lengths oq and 
02 as functions of m. 

We show in Fig. 2 oq as a function of m for different z. It is seen that for each z, there 
is a minimum value of Qq, around which there is a plateau that a^ changes slowly with m. 
We have pointed out in [10|, that for each z, there is also a plateau around the minimum 



of quark condensate varying slowly with m. The two m regions corresponding to the Qq 
plateau and quark condensate plateau almost coincide. 

And in Fig. 3, 02 as a function of ao for different z is shown in the (oq, 02) plane. It 
can be seen that in the (00,02) plane, there is a turning-point of the NJL 02(00) curve for 
each z, which corresponds to the minimum value of oq for each z. It is noticeable that all 
the turning-points for different z in the (09,02) plane are on a line! Remembering that Og 
and 02 are constrained to lie on a universal curve Eqs. ([121) with Lq = 0.58 ± 0.015 ni^^, we 
show this curve as double solid lines in the (oo, 02) plane. It can be seen that the turning- 
point of the NJL 02(00) curve for z = 0.7 falls on the universal curve. This turning-point 
(oo = 0.1796 m^^, 02 = —0.0489 m^^) is shown as a sohd circle in the (oo, 02) plane. 

To compare our results with other theoretical predictions, we also plot all the predictions 
(in units of m-^) in Fig. 3: 1), The Weinberg values (oo = 0.158, 02 = -0.045) H; 2), 
The chiral perturbation theory at one-loop (ChPT Hoop) |l5l predicted (oo = 0.200, 02 = 
—0.043); 3), The chiral perturbation theory at two-loop (ChPT 21oop) |16| predicted (oo = 



0.217,02 = —0.0413); 4), Olsson's dispersion-relation constraint of the heavy baryon chiral 
perturbation theory [jT^ resuhs (HBChPT Olsson) |g (ao = 0.235, as = -0.031). It can 
be seen that the turning-point of the NJL a2(cfo) curve for z = 0, i.e., (oq = 0.1576, 02 = 
—0.0446) shown as a x coincides with the Weinberg values shown as a star. 

V. CONCLUSIONS AND DISCUSSIONS 

Our results show that, in the framework of SU(2) NJL model, there is a turning-point of 
the 02(00) curve for each z in the {ao,a2) plane, and the turning-point corresponds to the 
region of the quark-condensate plateau; And the turning-point for z = is almost the same 
as the Weinberg values; Moreover, we find that all the turning-points in the (oq, 02) plane 
are on a line, which has a intersecting point with the universal curve at the turning-point of 
z = 0.7. 

We then can read the four parameters corresponding to this point: m^ = 8.22 MeV, 
A^ = 667.5 MeV, A^ = A/ x 0.7 MeV, and GAj = 4.36. And also we can calculate 
the other corresponding quantities at this intersecting point: the constituent quark mass 
m = 400 MeV and the quark condensate < qq >^/3= 223.4 MeV. 

Also we find that our results of pion-pion scattering lengths determined by the universal 
curve, i.e., the turning-point oi z = 0.7, are smaller than other predictions. The reasons 



maybe lie in that: a). From [13] we know that the Lq integrand of the universal curve include 
dominantly the p's contribution and other higher resonances. Here, we only consider the 
sigma exchange contributions. If one wants to investigate p's contribution, the extended 
NJL model should be used like in Ref. [|r^; b). We only calculate to the lowest order in 1/A^c 
expansion of pion-pion scattering Feynman diagrams. 
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FIG. 1. Feynman digrams for leading and subleading order quark self-energy niH and 6m, and 
meson polarization functions 11^^ and 5Il\f' . The solid and dashed lines indicate the quark 
internal meson propagators respectively. 
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FIG. 2. oo as a function of m for different z. 
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FIG. 3. The (00,02) plane include 1), a^ as a function of ao for different z i.e. 02(00) curve in 
NJL model; 2), the universal curve shown as solid double lines; 3), other theoretical predictions. 



